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Abstract 

Heavy Quark Effective Theory (HQET) computations of semi-leptonic decays, e.g. 
B — )■ nil', require the knowledge of the parameters in the effective theory for all com- 
ponents of the heavy-light flavor currents. So far non-perturbative matching conditions 
have been employed only for the time component of the axial current. Here we perform 
a check of matching conditions for the time component of the vector current and the 
spatial component of the axial vector current up to one-loop order of perturbation the- 
ory and to lowest order of the 1/m-expansion. We find that the proposed observables 
have small higher order terms in the 1/m-series and are thus excellent candidates for a 
non-perturbative matching procedure. 
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1 Introduction 



B meson decays are an excellent source of information for constraining physics beyond 
the Standard Model. Precision based on a solid theory and advanced experiments 
is becoming increasingly important as we know that effects due to fields which are not 
present in the Standard Model are small. Next to leptonic decays, exclusive semileptonic 
decays are easiest to treat in theory. Take for example the decay B — )• ttIu which is 
relevant for a determination of V^b- Theory only needs to predict two form factors (in 
practice a single one dominates) from non-perturbative QCD. This is a strong motivation 
to extend the HQET programme of the ALPHA-collaboration fTI|6| to include matrix 
elements of all components of the weak heavy-light currents. And it is a significant 
step beyond what has been achieved so far, where only the HQET action and the time- 
component of the axial current were determined non-perturbatively [2 , 6] . 

Instead of the previous five we now need 19 parameters in order to have the effective 
theory defined non-perturbatively including all 1/m terms, namely all terms of mass 
dimension five in the action and dimension four in the currents. Therefore, 19 matching 
conditions are needed. It is important to choose them well. Each matching condition 
simply consists of a matching observable $i which is evaluated in QCD and in HQET 
— in the latter theory including the terms of order 1/m and no more. Setting ^^"""^ = 
^HQET (jg^gj-j-Qjjjgg fact defines) the parameters in HQET. What does it mean to 
choose the matching observables well? Ideally we would like each one of them to be 
sensitive to a single parameter in HQET, in practice we would like them to receive little 
contributions from terms of order 1 /m? in the effective theory. If such contributions from 
0(l/m^) terms are unnaturally large, they affect the determined parameters and then 
inflict unnaturally large 1/m? terms into the observables that one wants to determine 
from HQET after the matching has been carried out. One thus better chooses the 
matching observables in QCD which are strongly dominated by the terms of order mP 
and m~^. Since the ALPHA strategy consists of matching in a finite volume with 
Schrodinger functional boundary conditions, the size of different terms in the expansion 
is given in terms of z~" = {Lm)~^ with L the linear extent of the finite volume. 

Of course, in the whole process, the most important terms are those which appear at 
order m? , the static terms. They are simply dominating numerically. It is thus of impor- 
tance to make sure that those matching observables which determine the normalization 
of the static currents are chosen well. Due to the breaking of relativistic invariance 
we need to normalize the space and time components of the currents separately. Thus 
we consider the axial vector current Aq,A]^ and the vector one Vq^V^. Previously, the 
normalization factor Z^^"^ of Aq has been studied in detail j2|[5 -j? 7 - 9j . It is defined 
through a Schrodinger functional two-point function [2]. Since in static approximation 
^0 and Vk are related through the spin symmetry (see Sect. [2] for a more precise state- 
ment), the natural condition for Zy^^"*" follows from a simple spin rotation. However, 

Z^^"^ and Zy^^"^ do not appear in the Schrodinger functional two-point functions 
which have been considered so far. We are thus lead to either consider two-point func- 
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tions with more complicated kinematics or three-point functions. 

In fact three-point functions appear naturally, since they are also used to determine 



the desired form factor for B — t- ttIu 10-16 . One thus uses a process in the finite volume 
matching which is related to one of the desired infinite volume matrix elements and there 
is even a potential that higher order in 1/m terms cancel between the matching and the 
physical matrix element. On the other hand, these functions have not been considered 
before. We therefore evaluate them first in perturbation theory, including the one-loop 
parts. We can then verify that they are indeed dominated by the first two terms in the 
1 /m-expansion. 

The perturbative study is rather straight forward, since one of us has developed 
"pastor", a tool to carry out one-loop computations of Schrodinger functional corre- 
lation functions in a largely automatic manner. Still, the scope of this paper is not to 
consider the full system of 19 unknowns, but to study the two numerically dominating 
matching conditions for Zy^^"*" and Z?^^"*". The pastor software package was first 



introduced in 17 and the publication of a more thorough description along with the 



source code is planned for the near future. 

2 The large mass limit of QCD: Heavy Quark Effective Theory 

We consider QCD with at least three flavors, one of them massive, mb = m, and the 
others massless, in particular m^^ = md = 0. A pseudo-scalar state with the flavor 
content bd is written | Pj^^j , L) , with L denoting a single external (kinematical) length 
scale. Analogously a light pseudo-scalar state is \P^^ , L) and vector states are labelled 
with V instead of P. We are interested in matrix elements 

A^QCD(L,m) = (X„a,L|J."^x)|Xba,L), (2.1) 

of the QCD heavy-light current operators which correspond to the classical field 

jf{x) = Zji^^{x)r,Mx). (2.2) 

In particular we consider the axial vector current, = Ai,, with = 'j^ju and the 
vector current, Ji, = with T^^ = lu- In physical processes, L is an inverse momentum 
scale, but we will later use states in a finite periodic Lx Lx L volume. For the moment 
the relevant point is that L is the only scale apart from m. Then there is a perturbative 
expansion 

M^^^{L,fn) = {M^^^Y^\z)+g\L){M^^^)^^\z) + 0{g\L)), z = Lfn . (2.3) 

We will specify the renormalization scheme for g, fn when it becomes relevant. The 
renormalization factors Zj of the flavor currents are to be chosen such that the currents 



satisfy the chiral Ward identities 118, 19 . In the large mass limit, m — )• oo, L fixed, the 



matrix elements Ai^^^ are logarithmically divergent 20 



21 



(7WQC°)«(z) -7olog(z)F(o), 70 = -l/(47r2), z = Lfn. (2.4) 
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This hmit of QCD is described by an effective field theory, HQET. Up to corrections 
of order l/z, it is the static effective theory ^22j where the b-field is replaced by a 
two-component static field, 

Mx) ^ M^) = IC^ + ^q)Mx) , (2.5) 

with Lagrangiar{^ 

^stM = M^KSm + Do)Mx) ■ (2.6) 

The mass counter term 6m does not play a role in the following. The static flavor 
currents are form-identical with the QCD ones, for example Vq^^^{x) = '0u(^)7o'0h(3;)5 
j4^*'^*(x) = V'u(^)757fcV'h(3;). Chiral Ward identities fix the relative normalization of 
the static vector and axial vector currents but not the overall normalization. Further- 
more space and time-components are to be treated separately and the currents have 
an anomalous dimension in the effective theory. Choosing the lattice regularization we 
can in a first step define finite currents by renormalizing them in the lattice minimal 
subtraction scheme. The renormalized currents are then 

(j£f ),(x;^) = Ziat(/ia,5o) Ju"\x) = Zi,t{f,a,go)i^^{x)T,Mx) , (2.7) 
with a renormalization constant 

^iat(/"a, 5o) = 1 - 70 log(a^)5o + ©(fo) ; (2.8) 

which is common to all currents (see [23| for a pedagogical introduction). Their matrix 
elements 

7Wt*(L,/i) = Ziat(^a,<7o)(^udl"/^*"*(x)|Xba)stat , (2.9) 

are then finite. When we set fi = fn, they are equal to the corresponding QCD matrix 
elements up to higher order terms in 1/m, 

A^Q^°(L,m) = CZ'^\-g\fn))MT\L,fn) + 0(l/m) , (2.10) 

and up to the finite renormalization factor 

CZ''\g'') = l + Bj^g^ + Oig^). (2.11) 

The one-loop coefficients are 

Bvo - Baq 
Bau - Bvo 

^We are in the frame where has 
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= -0.137(1), (2.12) 
= 0.0521(1) = Bv,-B A,, (2.13) 
= -0.016900. (2.14) 

spatial momentum zero and HQET at zero velocity applies. 



Here eq. (2.12), due to [81124), and eq. (2.13), due to 25 , depend on the lattice regu 



larisation. They are given for the Eichten-Hih lattice action for the static quark, the 
0(a)-improved Wilson action for the light quarks and the plaquette gauge action. We 



note that eq. (2.13) follows from requiring a chiral Ward identity. On the other hand the 



bare currents Vq and are related by the spin symmetry of the static effective theory 



which is exact in lattice regularization. The difference, eq. (2.14), is therefore known 
very precisely from continuum perturbation theory 26 . Of course the renormalization 



of the fields and therefore in particular Bj^ are independent of the states in eq. (2.1 ). 



3 Matching conditions 

3.1 Definitions of correlation functions 

As discussed in the introduction, in the ALPHA strategy we use finite volume ma- 
trix elements to define the matching of HQET and QCD. These matrix elements are 
constructed in the Schrodinger functional, where they are exactly related to ratios of 
correlation functions, see [27] for more details. Here we define those correlation functions 
and ratios which are suitable for the matching of Vq and A^. 

We choose the Schrodinger functional with vanishing background field, denote the 
time-extent by T and the space-extent by L. As a shorthand we introduce (non-local) 
boundary fields 



L3 



J]C.(x)ro(y), o^(r) = ^3 j;cKx)rcj(y), 



x,y 



(3.1) 



where the first one creates a meson with flavor content ij at time zero and the second 
annihilates a meson with flavor content ji at final time T. The boundary quark fields 
Ci , Ci are defined in 



28 



For simplicity and because more sophisticated choices seem 
unnecessary, we take each flavor to have the same periodicity phase 9 in the boundary 
conditions ^/'(x = Lk) = e^^ip{x) , ■iIj{x = Lk) = e~^'^ip{x) . 

With these preliminaries we define boundary-to-boundary correlation functions (re- 
member z = fnL) 



Kf{6) 



-Uo'^MOudiik)) , 



and three-point correlation functions with the desired currents 



Fvo(xo;6',2;) 
JI^{xo;6,z) 



(Odu(75)^0'^^x)Obd(75)), 

4(0:i,(7iX(x)Obd(75)). 



(3.2) 
(3.3) 
(3.4) 



(3.5) 
(3.6) 
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3.2 Possible matching observables for Vo,Ak 



The defined correlation functions are easily combined to form the desired finite volume 
matrix elements, 



QCD 



{L,m) 



-Zy 
-Za 



Fy,{T/2;9,z) 
Jl{T/2;e,z) 



(3.7) 
(3.8) 



where we set T = L. As explained in 27 these ratios are equal to the matrix elements 
eq. (2.1) with the finite volume states such as |-Pba, L), all normalized to unity. We 
here neglect 0(a)-improvement, but this is used in the perturbative computations in 
Sect. H 

We now have good candidates for matching conditions which we write in the form 



$Q^°(L, m) = (L, fn) + log {Cr'^' (^'(m)) } + 0(l/m) , 



(3.9) 



with ^^^^ = log [L'^M^^^] . In this way the log(C™'^*'^'^)-term appears additively, 



which is advantageous once the 1/m-terms are included [2]. 
3.3 Checking their quality 



Expanding eq. (|3.9|) in the coupling we have 

: (Cl>^f)(0)+O(l/Z), 

: {^%-^)W + Bj^-^olog{am) + Oil/z) 



($QCD)(0)(,) 



The one- loop part can be rewritten as in eq. (2.4), namely 



with 



H 
H 



(0) 

J. 

(1) 
J. 



($QCD)(i)(^) + 70 log(z) = + 0(1/ z) 
= {^l-^)W + Bj^-jolog{a/L), 



(3.10) 
(3.11) 



(3.12) 

(3.13) 
(3.14) 



where we subtract the logarithmic singularity in z from ($J^°)(i)(z) such that H^jJ 
represents the one-loop coefficient of the matched static matrix element at renormal- 
ization scale 1/L. In this form the size of 1/m terms is directly visible as deviations of 
the left hand side of eq. (3.12) from We want to investigate these deviations in 

the following in order to ensure that eq. (3.7) and eq. (3.8) are good observables for the 
matching. 
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4 One-loop computation 



All the required quantities {Fyg, J^^, F^'^, Kf^, and their static counterparts) were 
calculated at the one-loop level using the pastor software package for automated lat- 
tice perturbation theory calculations f[7\. As input, pastor accepts a rather general 
class of lattice actions and observables defined in the Schrodinger functional. It will 
then automatically generate computer programs for the evaluation of all contributions 
of the observables under investigation up to one-loop order including improvement- and 
counter-terms. We did implement full 0(a)-improvement, including the terms propor- 



tional to arriq not written in eq. (3.7) and eq. (3.8). 



For the quantities in QCD, we choose lattice resolutions of L/a up to 40, while 
for the HQET counterparts lower resolutions up to L/a = 30 are sufficient to obtain 
reliable continuum extrapolations, c.f. Fig. [T} To determine the continuum limits, we 
employ the method described in |29| using the implementation provided by pastor. We 
choose d e {0, 0.5, 1.0} and z e {4, 6, 8, 10}. 
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-0.03 - 
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-0.04 



-0.045 - 
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-I 



(1) 

Vo 



gV'^(z = io) 



0.002 0.004 0.006 

(a/i)2 



0.008 



0.01 



Figure 1: Continuum extrapolation of i/vo Gygat one-loop level, 9 = 0.5. The round-off 
errors on the data points at finite L/a and the uncertainty of the continuum extrapola- 
tion for the static point are much smaller than the symbol size. 



We employ the mass-independent lattice minimal subtraction scheme 28 in which 
the 0(a) improved renormalized mass at scale // = 1/L is given by 



m(L) = Zin,iat(5o'a/L)mq [l -ha6m(5o) 



rua 



nin 



mo 



rur 



(4.1) 
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in terms of the bare mass of the lattice theory. At one-loop order we have 30,31 



-0.5 - 0.07217(2) Ci.5o + O(<7o), 
-^\og{a/L)gl + 0{gt). 



(4.2) 
(4.3) 



All calculations in pastor are performed with z = m{L)L as input. It inverts eq. (|4.1 



to obtain mo = ttiq''^ + QQiri)^' + 0(go) ^-i^d calculates the series 



O 



m, 



(0) 



(1) 



0(0) I m, 



+ 5o 



(0) 



mo 



O(o) 



m, 



(0) 



+ 0{gt) (4.4) 



for a given observable O(mo). For the evaluation of the diagrams of a Schrodinger 
functional observable, it is beneficial to work in a time-momentum representation. Due 
to the periodic spatial boundary conditions one does not have to perform a momentum- 
integration but a sum over a finite set of allowed lattice momenta of size {L/a)"^ . The 
round-off errors introduced by the numerical evaluation of this sum are estimated from 
the difference of long double precision and double precision results for representative 
parameters. Apart from this test we use double precision since it is roughly a factor 
three faster. The execution time to evaluate the numerically most challenging loop 
diagram at L/o = 40 was about 50 hours on a single core CPU (Nehalem). 



5 Results 



5.1 Tree- level 

We start the discussion of our results with the tree-level functions Gy^^(z) = {^y^^)'^^\z) 
and G^^_^{z) = (<I>^^^)^'^^(z). Together with the static values -ffj^^ they are displayed in 
Fig. H and Fig. || for three different values of 9. Curves are fits of the form H^^\l + 
hi/ z + h2/ z"^), fitted to the data with weights w{z) = l/z^. The fits are thus dominated 
by the results at large z. The coefficients hi, listed for the different cases in Table [l} are 
small. For all considered values of the 1 /m-expansion is well behaved and we can also 
be confident that the fitted coefficients are close to the true Taylor coefficients. Obvi- 
ously, from the point of view of tree-level, one would prefer ^ = where G^j^ {z) = Hj^J 
holds exactly. 



5.2 One-loop 

We get more information at one-loop order. In order to have all finite pieces defined, 
we need to specify the renormalization scheme for the quark mass. As stated in Sect.|4j 
we take m to be the renormalized mass in the lattice minimal subtraction scheme at 
scale n = 1/L. The continuum limit is taken as described in the previous section. 
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Figure 2: G^^^ 
the symbol size. 



(*r)'°'(^) 



in the continuum hmit. Errors are much smaller than 



e 


0.0 


0.5 


1.0 




hi /i2 


hi /i2 


hi /i2 




0.00000 0.00000 
0.00000 0.00000 


0.77621 1.11933 
-2.30791 1.57043 


1.05083 1.53061 
-3.06017 3.11951 




hi fi 


hi fi 


hi fi 


'^Ak 


0.05245 -0.00132 
0.03099 0.01054 


0.12513 0.00139 
0.10547 0.01225 


0.21893 0.01547 
0.19391 0.02929 




0.15093 -0.00923 
0.12100 0.00692 


0.08803 -0.00811 
0.07042 0.00139 


0.04548 -0.01340 
0.05268 -0.01728 



Table 1: Fit coefficients for GA^and Gvo- The upper row of fit coefficients for the 
one-loop results comes from the fits omitting the data at z = 4. 



The combination G^j^{z) 



eq. (3.12), is shown in Fig. 4 and Fig. 



.0, 



fit to the one loop data employing a function of the form 

G'^j^{z) = H^^J + hi/z + filog{z)/z, 



We perform a 
(5.1) 



choosing in this case constant weights, as only few data point are available anyway. It 
is compared to a fit of the same form, omitting the data at z = 4. The fit parameters 
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Figure 3: G^-,' = 



QCD^ 



(z) in the continuum Hmit. Errors are much smaller than 



for the one- loop quantities in Table [T] are not expected to be accurate estimates for the 
corresponding asymptotic expansion. The accuracy of the fits and the smallness of the 
coefficient fi, however, may be taken as an indication that higher order terms in the 
l/z-expansion are not very important for the considered range in z. 

The size of H^^^ is relevant for us only as a consistency check: for all cases it 
is a little smaller than the expected magnitude 1 / (47r) for a perturbatively accessible 
quantity. The interesting question is the magnitude of 1/z-terms as well as curvature 
when Gj^ are considered a function of 1/z. 

We observe that also at one- loop order the 1/z-terms in <I>j^ remain small, but 
= is not preferred any more. A choice 9 = 0.5 appears a good compromise between 
tree-level and one-loop. Take for illustration g'^ = A and z > 10 as it is typical in the 
non-perturbative application Then we roughly have a few per-mille 1/z correction 
at tree-level and an ~ 3% correction at one-loop. This is very acceptable. We then have 
all rights to expect that the 1 / z'^ corrections, which are omitted when HQET is treated 
non-perturbatively (4 ,32, are negligible and indeed the curvatures seen in Fig. |4] and 
Fig. [5] are small. 
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0.0 I ■ I 6* = 0.5 I • I 61=1.0 I A I 




Figure 4: Gy^{z) in the continuum limit, compared to the static result. 



0.0 I ■ I 6* = 0.5 I • I 6* =1.0 I A I 



y< -0.05 




Figure 5: G^^^(z) in the continuum limit, compared to the static result. 
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6 Conclusions 



The proposed three-point functions appear very useful. They are seen to be strongly 
dominated by the lowest terms in the 1/z expansion. As a consequence, the three- 
point functions may well be applied to fix the remaining two unknowns, Zy^^"^ and 
^HQET^ in the static approximation non-perturbatively. We would recommend 9 = 0.5, 
but the one-loop study does not suggest this choice to be much superior to ^ = 
or 9 = 1. At order 1/m the full system determining the 19 parameters has to be 



considered. Three of these parameters come from the HQET action 33 , two from the 
temporal components of the vector and axial vector current respectively and the spatial 
components of the currents require the inclusion of further six parameters each ^34j-|36] . 
A study of this system in perturbation theory is presently being carried out by the 
ALPHA collaboration. 

We can also confirm that the new package pastor is very useful in studying such 
problems in perturbation theory. This goes beyond issues related to the regularization 
such as renormalization factors or improvement coefficients. In fact, all results presented 
here refer to the z-dependence in continuum perturbation theory, since we were able to 
reliably take the continuum limit a/L — t- 0. We have presented the results in the lattice 
minimal subtraction scheme for the quark mass. They can trivially be connected to the 
MS scheme by using [sT] m(L) = (1 + 0.12228252) x m^{l/ L) + 0{f). 
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